Magnetic fields above the surface of a superconductor with internal magnetism 
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The author presents a method for calculating the magnetic fields near a planar surface of a 
superconductor with a given intrinsic magnetization in the London limit. He computes solutions 
for various magnetic domain boundary configurations and derives relations between the spectral 
densities of the magnetization and the resulting field in the vacuum half space, which are useful 
if the magnetization can be considered as a statistical quantity and its features are too small to 
be resolved individually. The results are useful for analyzing and designing magnetic scanning 
experiments. Application to existing data from such experiments on Sr2Ru04 show that a domain 
wall would have been detectable, but the magnetic field of randomly oriented small domains and 
small defects may have been smaller than the experimental noise level. 
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I. INTRODUCTION 

Starting with the discovery of ferromagnetic order in 
the superconductors HoMogSgii^ and ErRh4B^ about 
three decades ago, there has been increasing interest in 
superconductors with an intrinsic magnetization. In the 
above cases and some of the RNi2B2C compounds (R = 
rare earth), the magnetization is due to localized mo- 
ments of the rare earth ions and coexists with supercon- 
ductivity in some temperature range. In other materials, 
the conduction electrons may not only superconduct, but 
also carry some magnetization. A prominent example is 
Sr2Ru04^, which is believed to have a complex, time re- 
versal symmetry breaking p-wave order parameter;^ so 
that the orbital angular momentum of the Cooper pairs 
creates a magnetic moment. This is theoretically ex- 
pected to cause edge currents at sample boundaries and 
domain wallsi^ A similar effect has been suggested to 
occur in the d-wave superconductor Naa;Co02 • J/II2O.— 

The traditional experimental techniques for study- 
ing magnetic ordering phenomena are bulk probes 
such as muon spin rotation (/iSR) or (spin polarized) 
small angle neutron scattering (SANS). An alterna- 
tive approach is to use magnetic scanning techniques 
such scanning Hall probe microscopy (SHPM))^ scan- 
ning SQUID microscopyjiS magnetic force microscopy 
(MFM)iiii^ and magneto-optical techniques'^ These 
techniques measure the magnetic field some small dis- 
tance above the surface of the sample as a function of 
position. In many cases, their resolution, which is lim- 
ited by the probe size and probe-sample distance, does 
not quite reach the length scales typical for the magnetic 
structure. For example, one finds evidence for oscillatory 
magnetic order with a sub-penetration depth length scale 
in HoMoeSgii^ and ErRh4B4i^ For larger wave lengths, 
the Meissner effect in the coexisting state would screen 
the field and thus suppress magnetic interactions and 
destabilize the magnetic order. A similar situation has 
been observed in the superconductor ErNi2B2C using 
SHPMJ^ The direct observation of the magnetic fields 
generated by edge currents or domain walls in Sr2Ru04 



is an ongoing effort which has not produced any evidence 
so fari^ii^. 



To plan such scanning experiments and to interpret 
the resulting data, it is important to understand to what 
extent the magnetic field generated by a spatially vary- 
ing magnetization inside the sample is propagated to the 
probe. In normal materials, this is a straightforward 
magnetostatics problem. In superconducting samples, 
however, it is complicated by the Meissner screening. 
In this work, I present solutions for an infinite planar 
sample surface by incorporating the presence of a mag- 
netization M(x) into a London model and solving the 
resulting equations, using a 2D Fourier transform. This 
approach follows earlier theoretical work on supercon- 
ductors with an internal magnetization due to localized 
magnetic moments fiiii^iiSii^ but should also apply if the 
magnetism is of different origin, such as a spin or or- 
bital moment of the Cooper pairs. To consider the prac- 
tically very likely case where the limited measurement 
resolution leads to significant averaging over several do- 
mains or other features with some degree of randomness, 
I present a spectral analysis. The resulting relations were 
employed for the analysis of the data in Ref. [ij and may 
also be used to analyze null results where no field vari- 
ation is detectable at the experimental noise level. A 
similar analysis may also be useful for interpreting sur- 
face sensitive ^SR experiments, which only average over 
a thin layer at the sample surface^!. 



The paper is organized as follows: In section [TTl I de- 
rive the main equations of our model from a generalized 
Ginzburg-Landau (GL) functional. Those equations are 
solved in a general framework in section IIIII In section 
IIVI I discuss simple domain wall and dipole configura- 
tions as examples. Relations between the spectral densi- 
ties of the magnetization and the magnetic field in vac- 
uum are computed in section |Vl Section IVII applies the 
results to recent magnetic scanning worki^ on Sr2Ru04. 
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II. MODEL 

In order to fully describe the interplay between mag- 
netism and superconductivity in a phenomenological ap- 
proach, the magnetization M and the superconducting 



order parameter ip have to be computed self consistently, 
taking mutual interactions into account. This has been 
done in Refs. [T^lIllillS using the generalized GL func- 
tional 
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where B = V x A. In this work, I will assume this task 
to be partially solved by starting with a given M(x) and 
computing the resulting magnetic field, taking shielding 
currents into account. This approach is clearly justified if 
the magnetic energy scale is much larger than the super- 
conducting one, so that the effect of superconductivity 
on magnetism can be neglected. However, it is also rea- 
sonable if the result of a self consistent calculation for 
M is (approximately) known, for example from bulk cal- 
culations, and one is mainly concerned with the effect 
of the reduced screening at the surface on the observ- 
able field. The errors introduced by this treatment will 
then primarily be due to the effect of the surface and the 
modified screening on M. 

Writing the order parameter as V' = IV'Ig*'^ and intro- 
ducing the London penetration depth defined by — 
47r(|?/;p/m*)e*^/c^, variation of Eq. ^ with respect to 
A leads to 

V X (V X A - 47rM) + l/A^A = ($o/27rA2)V0. 

By performing a line integral over a closed loop after mul- 
tiplying by and using the Stokes theorem, one obtains 

V X (A^V X B) + B = 47rV x (A^V x M) + $of • (2) 

f is a sum of 2D 5 - functions representing vortex cores, 
which will be ignored in the following. This result can 
be obtained directly by treating the supercurrent density 
js as macroscopic current in the macroscopic Maxwell 
equation and thus substituting V x (B — 47rM) = 47r/c js 
into the London equation V x (47rA^/c js) -I- B = <I>of. 
Although the above is valid for a spatially varying su- 
perfluid density Ug = \ip\'^/m*, I will assume A to be 
constant, which leads to the more familiar form 

V X V X B + l/A^B = 47rV X V X M. (3) 

I assume the magnetic superconductor to occupy the 
lower half space z < 0. In vacuum [z > 0), the mag- 
netic field must satisfy V • B = and V x B = 0. At 
the interface, the normal component of B and the tan- 
gential component of H ~ B — 47rM must be continuous. 



Note that M enters Eq. [3] only through the microscopic 
current density jjyi = cV x M. Thus, one may also start 
directly from an intrinsic current density jM rather than 
M, which is more natural if an edge current is known 
from microscopic calculations, for example. The appear- 
ance of M in the tangential boundary condition can also 
be eliminated by replacing it with a discontinuity in M 
just below the surface. 



III. SOLUTION 

In Ref. [13, the field geometry of a vortex penetrating 
the surface of an anisotropic superconductor for a gen- 
eral orientation of the vortex and the main axis of the 
effective mass tensor with respect to the interface has 
been computed. In the vortex problem, the right hand 
side appearing in the London equation Eq. ([3]) is a 2D 
delta function instead of the magnetization term. I use 
the same technique, but only present the calculations for 
the isotropic case for the sake of simplicity. 

The Maxwell equations in vacuum can be satisfied by 
writing the magnetic field as B = — with — = 0. 
A suitable solution has to be matched to a solution of 
Eq. ([3]) at z = 0. I solve this problem using a 2D Fourier 
transform (FT) in the xy-plane, i.e. by writing a func- 
tion A{x,y,z) as A{r\\,z) = (27r)"2/ d2ki(k, z)e*''""ii , 
with r|| = {x,y) and k = {kx,ky). (Note that I use the 
same symbol for 2 and 3 dimensional vectors, implying 
that the z- component vanishes for the latter.) For z < 0, 
the field B is decomposed as Bq -I- Bi. Bq is a partic- 
ular solution of the inhomogeneous London equation ([3]) 
in full space with boundary conditions at infinity and 
with the right hand side suitably extended to z > 0. Bi 
is a general homogeneous solution chosen to satisfy the 
matching condition at the interface. Under the 2D FT, 
Eq. §1 transforms into {P + l/A^ - d^/dz^)Bi = 0, so 
that Bi(k,z) = BkC^^ with K = y^W + T/>^. In vac- 
uum, — = has the solutions <I>(k, z) = $ke~'^" so 
that at z > 0, B(k, z) = {-ik + kez)^iie^''^ . 



3 



Hence, V • Bi = together with the e^, k and k x 
components of the continuity conditions for B and H = 
B - AttM at z = read 



= ik • Bk + Ke^ ■ Bk 
A:$k = ■ fBk + Bo(k,0) 



-ifc$k = k- (^Bk + Bo(k,0) -47rM(k,0) 

= (kxe^) • (Bk + Bo(k,0)-47rM(k,0) 

The last equation for the in-plane transverse compo- 
nent of Bk is already decoupled and can be dropped if 
only the vacuum field is to be computed. Solving the first 
three equations for $k leads to 

k{K+k)<i>k = Ke,-Bo{k, 0)+i k- (Bo(k, 0) - 47rM(k, 0)) 

If the 2D FT of the inhomogeneous solution, Bq, can- 
not be obtained directly, the 3D FT Bo(q) of Bo(r) can 
be obtained from M(q) by solving the 3x3 linear system 

-q X (q X Bo(q)) + l/A^Bolq) = -4^q x (q x M(q)) 

The solenoidal condition V • Bg = will always hold as 
q • B(q) = A2 q • [q X (q X B(q)) - 4^q x (q x M(q))] - 
0. The 2D FT at z = is then simply Bo(k,0) = 
(l/27r) J dqz'BQ{q) with q = k + g^e^. While this ap- 
proach to the inhomogeneous problem is very convenient 
for numerical evaluation and can be generalized to the 
anisotropic case (cf . Ref . [23 ) , it is useful to derive an ex- 
plicit solution. The component of M parallel to q does 
not contribute to q x M. For components of Bg and 
M orthogonal to q, q • Bq — automatically and the 
vector products simplify to scalar multiplication. By de- 
composing M and Bq into components along the unit 
vectors e^y = q x Sx/jq x e^j orthogonal to q in the 
xy-plane and ej^ = q x e^y, the inhomogeneous solution 
thus simplifies to 



Bo.xy 



1/A2 + fc2 + 



^2 ^ L2 ^ ^a^T^^-y 



1/A2 + fc2 + g2 



(5) 



To evaluate Eq. ([l]), the inverse z-FT must be carried 
out in order to obtain the values at z = 0, unless M is z- 
indcpendent. It turns out that Bo.xy does not enter Eq. 
(m because of the dot product with k. Projecting i?o,i 
onto the z and k direction, substituting Eq. ([5]), and 
expressing everything in terms of Mz and k • M using 
M I = + ,1' „k • M leads to 



k{K + fc)$k = 2 dq. 



-2i dqz 



k^K-iqz) 
l/X^ + k^+ql 
K{K - iqz) 



l/A2 + fc2-f g; 



Af.(q) 
k.M(q) (6) 



The iqz terms in the numerator of the fractions can be 
dropped if one uses the convention that M is extended 
to z > as an even function so that those terms do not 
contribute to the integrals. Eq. ^ can be summarized 
qualitatively as follows: For in-plane components of M, 
the source of B outside the superconductor is the diver- 
gence of M averaged over one penetration depth below 
the surface. For the normal component, an additional 
derivative is taken, thus increasing the multipole order 
of the vacuum field by one. The small k components 
of the field are just those resulting from the subsurface 
magnetization and its image obtained by reflection about 
a plane A below the surface. 

I would like to point out that for the solution method 
to work as described, the interface must be planar and 
A2 should be constant. M on the contrary can be an ar- 
bitrary function. However, if V x M = so that Bq = 
solves Eq. ([2]), the requirement that A2 may not depend 
on z can be dropped at the expense of solving a more 
complicated ordinary differential equation instead of the 
Laplace equation to obtain Bi(k, z). For example, dis- 
continuities in A2 as a function of z could be treated 
by matching additional continuity conditions. A depen- 
dence of A2 on a; or y on the other hand would mix dif- 
ferent k components and thus generally forbid a simple 
analytic solution. 



IV. EXAMPLES 

A. Discussion of Table [I] 

As examples for various simple, representative configu- 
rations in M, I have calculated the field of ferromagnetic 
domain walls, where M changes sign, and dipoles at the 
surface. For z 3> A, simple approximate expressions in 
real space can be obtained. The results are shown in ta- 
ble [D The approximations are based on the fact that for 
z; 3> A, only fcx <C 1/A contributes to the Fourier integral 
Bz{x,z). To second order, 1/{K + \k\) w Ae"!''^!^ and 
K w (1/A)(1 -|- (fcA)2/2). Thus, the approximations in 
table U are good to three and two orders beyond lead- 
ing order in k for cases (1),(3) and (2), (4), respectively. 
In the following, I discuss the far fields obtained from 
these approximations. To understand those, it is useful 
to recall that for magnetostatic problems in the absence 
of macroscopic or supercurrents, — V2<I> = — 47rV • M. 
Thus, V • M acts as a magnetic charge by analogy with 
electrostatics. The field of a discontinuity of the in-plane 
component of M [case (1)] is just twice that of a mag- 
netically charged line with linear charge density Air 2M A 
situated A below the surface. It can be imderstood as the 
charge density due to the discontinuity of 2M in the mag- 
netization, which is screened by supercurrents over one 
penetration depth. The additional factor two formally 
comes from the extension of M to z > 0. For a disconti- 
nuity in the normal component Mz [case (2)], one obtains 
the dipolc field of two lines of opposite magnetic charge 
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with a charge density MA as above and a separation of 
A. Again, this can be understood as the screened field of 
the discontinuity in the magnetization occurring at the 
surface. 

The locahzed dipole pointing into or out of the surface 
[case (3)] has a a quadrupole field. An in-plane moment 
[case (4)] on the contrary has a dipole field to leading 
order. For a dipole chain, i.e. M = •fhS{x)S{y)ez, I 
obtain the same results as for a single dipole oriented 
in the z direction [case (3)] with m = Am. The same 
analogy can be drawn for case (4). 

It is also of interest to consider a configuration where 
a magnetization Me^ is localized over a width w around 
X = 0. This situation may be encountered at an antiferro- 
magnetic domain wall, where canting of antiferromagnet- 
ically ordered in-plane moments produces a local net out- 
of-plane magnetization. The corresponding exact solu- 
tion is a superposition of two solutions for normal discon- 
tinuities with opposite signs, shifted by w. If is smaller 
than all other length scales, M{x) can be replaced by a 
delta function: M{x) « {J M{x') dx')5{x) = wM5{x). 
The resulting field is then simply the derivative of that 
of a discontinuity in A/^, i.e. for \, 



Bz{x, z) 



(z + A) 



dx'^ (z + A)^ + 



B. Periodic configurations 

If the magnetization is periodic, the Fourier integrals 
turn into Fourier sums. If the period L is large {L > 27rz), 
features with size comparable to L can be resolved in each 
unit cell individually and look similar to a solution ob- 
tained from a constant continuation of M outside that 
cell. For shorter periods however, the superposition of 
many such single cell solutions largely cancels out. For- 
mally, this follows from the fact that the wave vector 
k takes only integer multiples of 2'k/L. Therefore, the 
dominant contributions at fc < \/z considered in table U 
are not present an the leading term becomes that of the 
lowest wave vector k = 2tt / L. This results in an expo- 
nential suppression by e"^'^^/^ of the lowest harmonic of 
the field variation at a height z above the surface and all 
higher harmonics being negligible for z > L. 



C. Effect of smoothing 

As the exponential cutoff in the inverse FT becomes 
increasingly sharper for larger z, the far field from any 
feature of finite size in the magnetization will always be 
determined by the lowest non-vanishing power of k in 
k^k for sufficiently large z. For a smooth domain wall, 
M(q) has less weight at large q compared to a sharp dis- 
continuity of equal magnitude, but the values at small q 
are affected little. Thus, the asymptotic results for sharp 



discontinuities in M are still valid if M changes smoothly 
over a width w as long as z ^ w. For example, the far 
field of a single domain boundary does not depend on the 
length scale over which the magnetization changes in the 
xy-plane, but only on the difference between the asymp- 
totic values of M on both sides. Hence, the approximate 
results in table U are of rather general validity. For a 
graphic illustration, see Fig. [21 which is discussed in Sec. 
fVll The situation is slightly different for a nontrivial z- 
dependence of M: the g^-cutoff is always determined by 
the 1/(1 -I- (Agz)^) factor for small k. In real space, this 
corresponds to an exponentially weighted average over a 
distance of A below the surface. 



V. SPECTRAL ANALYSIS 

The resolution of currently available magnetic imaging 
techniques is often not sufficient to resolve an actual do- 
main structure. In practice, averaging occurs both due 
to the imaging height z and the finite sensor size. Here, 
I will only consider the more universal height effect. For 
a height larger than the typical domain size, the mag- 
netic field represents an average over several domains. As 
shown in the previous section, this would lead to an ex- 
ponential suppression for periodic configurations. How- 
ever, domains will usually not be strictly periodic, but 
have some distribution of size. Therefore, a statistical 
description is most adequate. Assuming a given corre- 
lation function and thus spectrum of the magnetization 
M(r), I compute the spectrum of the resulting magnetic 
field. 

The spectral function of two functions /(r) and g{r) in 
d dimensions, Sfg{q) = / dV(/(r) ^(r r'))e-*'''i, sat- 
isfies the relation (/(q)*g(q')> = (27r)''(5'*(q - q')5/g(q) 
Here, (•) stands for an ensemble average over different 
realizations of / and g. It is implicit to this definition 
that the correlator (/(r) g{r + r')) is independent of r. 

It is convenient to introduce the propagation coeffi- 
cients 



7q 



2k^X'^-fq 

—2ikaKX^jq {a^x,y) 

(K - iqz)e-^'= 
A2(fc + A0(l/A2 + /c2 + g2) 



and rewrite Eq. ([S]), using i?z(k, z) = fc$ke '^^ 
Bz(zo,k) = / Cq,QMQ(q). It follows that 



(7) 
(8) 

as 



27r / dq^^C* ,,Cq,/35M,a,3(q) (9) 
27r / (jg^ |Cq,a|^5'A/^aa(q) 

47r j dq:, ^ Re(c* „Cq,^S'A/,Q/3(q))(10) 
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Case 


M(r) 


Bz (exact) 


Bz (approximate) 


(0) 


M{X) By 








(1) 


Msgn{x)ex 


4A/ f dk ifcxg-|fc|z 
J \k\ + K 


4M /dA; Ae"'"e-I'=l("+^) - 8M\, ^ to'*' o 
/ (z + \r + 


(2) 


Msgn{x)ez 




(l/A2+fc2)(jfc|+A-) 


4M /"dfc ( ifc)A2e^*=^e-l^-|(^+^) - SA/A^ , ^ + ^ , 
J ax (2 + A)^ + 


(3) 




2m 


f°° k 


2m r dk fc^AJo(fc.)e-'-^) = 2mA6(^ + " ^'^'Jf 


(4) 


me^(5^(r) 


2m f dk ^ ^ cosyjJi(fcr)e '^^ 
Jo k + K 


2mcos<z) A dk k^ Ji{kr)e~'^^-^^^^ = 2m. cos ip , „ '^'^/^ Al.r ;o 
A, "^(r^ + (z + A2))5/2 



TABLE I: Exact and approximate expressions for the magnetic field above the superconductor for domain walls and dipoles 
with an in- plane and out-of-plane magnetization. In case (4), ip is the angle between the x-axis and r. The approximate 
expressions in the last column are valid for z 3> A. 



I have used SM,af3 as a short hand notation for the spec- 
tral function SmuMP of two different components of M. 
Similar expressions can be written down for spectral den- 
sities involving other components of B. 

Because of the presence of the surface, the assumption 
of translational invariance in the z-direction implicit to 
the definition of Sm is by no means trivial. If the presence 
of the surface does not affect the structure of M too 
much, and the range of the surface influence is much 
shorter than A, most of it should average out because Bz 
is sensitive to what happens within a layer of thickness 
A below the surface. However, the z-invariance is only 
required in order to define spectral functions Sm.oP- If 
the interface (or other effects) do break the z-invariance 
of {Ma{r) Mp{r+r')) so that it depends on both z and z', 
it is still possible to derive an expression for Sb^ similar 
to Eq. (fTO|) . however involving a double integral over q^. 

A statistical analysis will be most relevant when the 
measurement height z is much larger than any of the 
intrinsic length scales of the variation of M, i.e. too large 
to resolve individual features. In this case, Sm{)<- + Qz^z) 
will not have a strong k dependence in the small k region 
surviving the e~'^^ cutoff and can be approximated by 
Sm{Qz^z)- For z ^ A, a similar approximation can be 
made for the propagation coefficients Cq^a and 7q: 

Cq,Q ~ -2ika\-fq {a = x,y) (11) 

-2fe(z+A) 

i7,r - (12) 

In any case, the properties of M only enter via the inte- 
grals 

/^9.^_^^.2;,2+g2;,2^M,.Mq) (13) 

As argued above, it will often be a good approximation to 
set k = 0. In many cases, Sm{<1z^z) will have a peak at 



some wave vector (and consequently at —qo), similar 
to the illustration in Fig. [TJ For simplicity, I assume 
that there is only one such maximum. A finite q^ is a 
signature of an oscillatory behavior of M. The width of 
the peak corresponds to the inverse coherence length of 
the oscillation or the correlation length for q^ = 0. 

The integral \V6\ can be approximated further in two 
limiting cases. If the coherence length of M along the z 
direction is much larger than A, then Sm{<1z^z) is sharply 
peaked, and the kernel 1/(1 -f g^A^) can be replaced by 
1/(1 + 90^^) and pulled out of the integral. In the op- 
posite limit, the peak in Sm{(1z^z) is much wider than 
1/A so that the qz dependence of Sm can be neglected 
entirely and one obtains 

( d^z-—^-YT^SM.ai}{q) ~ (7r/A)S'A/(0) 

J 1 + qi^ 

Assuming that the a-component of M dominates, the 
respective diagonal term of Eq. (flOl) takes the form 

for a = x,y. If the Mz component is dominant, (fc^A)^ 
in the prefactor must be replaced by (/cA)^. Similar ex- 
pressions can be written down for the off-diagonal com- 
ponents of SM^af}- For an illustration of the relation be- 
tween a measured B, and its spectral function, the reader 
is referred to Ref. [ij- Integration leads to simple expres- 
sions for (B^) = (27r)-2/d2k5'5jk): 

(Bl) = T(7TA)^^*^^""^'^) fora-x,y (14) 

(B'z) = ^(^^^M.-(O)- (15) 

Those expressions can be used to estimate the sig- 
nal expected in a scanning experiment or to estimate 
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FIG. 1: (Color online) Illustration of different one- 
dimensional domain structures for |M| = 1. (a) Correlator of 
M, and small sample of the corresponding real space domain 
structure (inset), (b) Corresponding spectral functions. The 
domain width in the red (lower) structure is exponentially dis- 
tributed with mean 1, which gives {AI{x)M(x + x')) = e^^'^ 
and SM{k) = 1/(1 + fc^/4). The domain width in the blue 
(upper) structure has a Gaussian distribution with mean 1 
and standard deviation 0.3, leading to an oscillatory correla- 
tor and peaks in the spectral function at fc = ±7r. The curves 
for this case were obtained numerically from an ensemble of 
10^ domains. The dashed lines in (b) represent the propaga- 
tion factors e^'*^'"^ for heights z = 2 and z = 0.2. For large z, 
only SM{k ~ 0) is relevant. 



Sm{Q) from the observed field variation. Note that if 
the variation of the domain size is sufficiently large for 
(M(r)M(r + r')) to be essentially non-negative, as for 
the exponential width distribution in Fig. [Jl SM,aaiO) = 
/ d^r' {Ma{r)Ma:{r + r')) can be interpreted as the 
product of a correlation volume / cPr' {Ma{r)Ma(r + 
r'))/(M^) and the mean square magnetization (M^). 



VI. APPLICATION TO Sr2Ru04 

Based on various evidence, it is believed that Sr2Ru04 
is a spin triplet superconductor with a p-wave order 
parameter of the same symmetry class as ± ikyi^^ 



Convincing evidence that the order parameter is indeed 
time-reversal symmetry breaking (TRSB) has recently 
been obtained by Sagnac-interferometry experimentsi^ 
Such a TRSB order parameter is expected to cause chiral 
currents at sample edges, domain walls or impurities i^i^ 
The direct observation of such effects in Sr2Ru04 is an 
ongoing effort. So far, the most direct indication of 
spontaneous fields is given by /iSR data, reporting "a 
broad distribution of fields arising from a dilute dis- 
tribution of sources" Phase sensitive tunneling mea- 
surements support the notion of small chiral domainsi^ 
Scanning Hall probe and scanning SQUID microscopy 
experimentsi^ii^'^l on the other hand did not detect any 
sign of a spontaneous magnetization associated with su- 
perconductivity. 

The magnetic scans of the a6-face in Ref. [H showed 
neither localized features nor a random field variation 
that could be attributed to TRSB. Even holes that were 
drilled using a focused ion beam (FIB) failed to show a 
magnetic signature. Lacking suitable theoretical models, 
a quantitative analysis of those null results was inconclu- 
sive. In this section, I will use the results derived above 
to compute the expected field from domains and defects. 
This allows to set certain limits on the internal magneti- 
zation strength that would be consistent with the data. 

A complete description of chiral domains requires a 
self consistent computation of the order parameter and 
magnetic fields. This has been carried out using micro- 
scopic theory^ and a GL approach^ without considering 
the effect of a surface. Such detailed calculations gen- 
erally require a numerical solution. Taking the presence 
of a surface into account leads to a further complication. 
Thus, they are rather cumbersome for the purpose of 
data analysis and planning experiments. 

For a domain wall along the a; = plane in an infinite 
sample, the results of those computations generally show 
a current along the y-direction which decays over about 
one in-plane coherence length ^ab in the a;-direction and 
changes sign. The counterflowing current decays on the 
scale of A such that the magnetic field far inside each 
domain vanishes. Indeed, one can obtain good fits of the 
form 

= 73|/X^sgn(^)(^"'-''^' - e"'"'^^^) (16) 

to the numerical results for the magnetic field of Ref. @, 
with A = 2.2 ^ab, I = 1.5 S,ab, and Bq = 87 G. I used ^ab = 
66 nm and A = 150 nni^ to compute the thermodynamic 
critical field entering the prefactor of the result of Ref. 
0. It is easy to show that this expression is a solution 
to the London equations for A = A in the presence of a 
chiral current density jy{x) = — (ci?o/47r^)e~l^l/^. This 
current can be identified with an internal magnetization 
Mz = (i3o/47r)sgn(x)(l - e^l"'!/'^) in the z direction due 
to the orbital magnetic moments of the Cooper pairs. 
The value of A is in good agreement with the value of 
K = \/iab = 2.5 assumed in Ref. 0. This suggests that 
the London approach captures screening effects quite ac- 
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curately, and therefore should give a good approxima- 
tion for the field above a surface. Of course, this will 
neglect some features in the full GL solutions. For exam- 
ple, the latter show a slight depression of the superfluid 
density near the domain wall, which is neglected here by 
assuming a constant A. Similar modifications of both the 
superfluid density and the chiral current density at the 
surface should be small due to the short c-axis coherence 
length = 0.05 £_ab- The anisotropy is of no consequence 
because all currents flow along the a6-plane so that only 
Xab matters. 

Since the far fleld of the domain wall only depends on 
the difference of the asymptotic values of AI away from 
the domain wall, and experimentally ^ab — 66 nm < A = 
150 nm <C 2 ~ 1 ^m, it is appropriate to use the large- 
z results for a sharp discontinuity in M to analyze the 
results of Ref. [H. If accurate results at a z ^ 3A are 
required, Eq. ([4]) together with the above approximation 
for Bo^z{x) should be used. Fig. [2]shows the field profiles 
for those two methods and the exact result for a sharp 
domain wall at different heights. 
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FIG. 2: (Color online) Field profile at diff'erent heights z 
above a domain wall were changes sign. The curves la- 
beled "smooth" and "sharp domain wall" were obtained from 



the full solution with M^{x) = Msgn(a:)(l — e 



) and 



Mz{x) = Msgn(a:) respectively. The latter solution and the 
corresponding large z approximation as shown are given by 
table m case (2) . The value of the domain boundary width 
was chosen according to the fit to the numerical results for 
Sr2Ru04 from Ref. [y. z = 8X corresponds to the Hall probe 
scans from Ref. [la. 



In the experiments, the rms noise levels were 35 mG 
and 0.45 mG at imaging heights z = 1.2 //m and 2 
/im in the Hall probe and SQUID scans, respectively. 
Correction factors due to oversampling and averaging 
over the size of the SQUID pickup loop are of order 
unity. An isolated domain wall would result in a field 
of 5 • 10-347rM = 0.4 G and 2 • IQ-^AnM = 0.17 G re- 
spectively at the above imaging heights and Bq = 87 G, 
and should be clearly visible in the data. Thus, there 
was either no such domain wall in the scanned area, or 
its magnetization was 47rM < 7 G and 47rM < 0.2 G 
respectively, so that it was hidden by sensor noise. 

The exact calculation of the signature of a hole is 
more difficult because the translational invariance of the 
boundary conditions is broken. However, if the diameter 
of the hole or defect, which I assume to extend along the 
z-direction normal to the surface, is much smaller than 
A, the absence of superfluid in it can be neglected and the 
dipole calculation should be a good approximation. For 
a hole or defect with a volume of ^j^^, the maximum field 
according to case (3) in table [J at z = 1.2 fim and 2 /im 
is 1.2 • IQ-MttM = 1.1 mG and 2 • IQ-^AnM = 0.17 mG, 
respectively. This signal would be nearly undetectable 
at the experimental noise level of Ref. [l^. Furthermore, 
the extent of the defect along the z-direction could be as 
small as — 0.05 ^ab, and the order parameter is not 
necessarily suppressed entirely. For a columnar defect on 
the other hand, one factor of ^ab has to be replaced by A, 
decreasing the limit on M only by about a factor three. 
Since the FIB drilled holes in the experiment were sig- 
nificantly larger (about 1 /im), they have both a larger 
moment and less effective Meissner screening. This leads 
to a stronger signal whose calculation goes beyond the 
scope of this paper. 

One can also estimate the signal expected from a ran- 
dom configuration of small domains. The smallest con- 
ceivable domain volume is on the order of ^^i,£,c- Assum- 
ing that the domain size fluctuates enough to use this as 
correlation volume, Eq. p5p implies that the rms signal 
could be as small as 0.7 mG and 0.16 mG, again less 
than the experimental noise. A domain size distribution 
that does not satisfy the assumptions leading to Eq. ([15]) 
may result in even smaller signals. Thus, the possibility 
of very small domains cannot be ruled out. 

In all the above cases, the smaller imaging height of 
the Hall probe does not compensate for its large noise 
compared to the SQUID. Assuming the predicted mag- 
nitude of the chiral currents^ is correct, the calculations 
show that any domain wall should have been detected 
by the measurements. Small defects on the other hand 
might easily have been hidden in the noise. It is also 
possible that a random signal from domains would have 
been too small to observe, especially if the domains are 
short in the c-direction or very homogeneous in size while 
not too large. However, it appears that one should not 
take the notion of a magnetization due to p-wave pairing 
too literally. It was shown^^ that the chiral currents can 
in general not be written as the curl of a global magncti- 
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zation. and that the effective value of M depends on the 
type of domain waU. 

Note that a similar analysis of the magnetic scanning 
data of Refs. [TsIITgI . also making use of the relations 
derived in Sec. Illli has been carried out in parallel with 
the present work in Ref. [l^. 

VII. CONCLUSION 

I have presented a model for a superconductor with 
an intrinsic magnetization by combining the macro- 
scopic magnetostatic Maxwell equation with the Lon- 
don relation and obtained the field geometry at a pla- 
nar superconductor-vacuum interface for a given spatial 
variation of the magnetization. Solutions for a range of 
specific magnetic domain boundary configurations give 
simple expressions in the limit of a large height above 
the sample. If the height above the surface at which the 
magnetic field can be measured exceeds the characteristic 
length scale of variations in the magnetization, a spectral 



analysis can be used to relate the spectral densities of the 
two at resolvable wave vectors. If a specific model for the 
structure of the magnetization is at hand, a direct com- 
parison with the measured field is possible. Otherwise, 
some simplifying assumptions give a simple estimate re- 
lating the spectral density of M at the superconductor - 
vacuum interface to the observable spatial rms-variation 
of the magnetic field. As an example for an application, 
I have applied my calculations to recent experimental re- 
sults on Sr2RuO,f^, concluding that large chiral domains 
would have been visible in those experiments, but small 
domains and defects may have escaped detection. 
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